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Abstract: It has been known for a long time that vacuum polarization in QED 
leads to a superluminal low-frequency phase velocity for light propagating in curved 
spacetime. Assuming the validity of the Kramers-Kronig dispersion relation, this 
would imply a superluminal wavefront velocity and the violation of causality. Here, 
we calculate for the first time the full frequency dependence of the refractive index 
using world-line sigma model techniques together with the Penrose plane wave limit of 
spacetime in the neighbourhood of a null geodesic. We find that the high-frequency 
limit of the phase velocity {i.e. the wavefront velocity) is always equal to c and 
causality is assured. However, the Kramers-Kronig dispersion relation is violated due 
to a non-analyticity of the refractive index in the upper-half complex plane, whose 
origin may be traced to the generic focusing property of null geodesic congruences 
and the existence of conjugate points. This puts into question the issue of micro- 
causality, i.e. the vanishing of commutators of field operators at spacelike separated 
points, in local quantum field theory in curved spacetime. 



1. Introduction 



Quantum field theory in curved spacetime is by now a well-understood subject. How- 
ever, there remain a number of intriguing puzzles which hint at deeper conceptual 
implications for quantum gravity itself. The best known is of course Hawking radia- 
tion and the issue of entropy and holography in quantum black hole physics. A less 
well-known effect is the discovery by Drummond and Hathrell [1] that vacuum po- 
larization in QED can induce a superluminal phase velocity for photons propagating 
in a non-dynamical, curved spacetime. The essential idea is illustrated in Figure 1. 
Due to vacuum polarization, the photon may be pictured as an electron-positron pair, 
characterized by a length scale Ac = m~^, the Compton wavelength of the electron. 
When the curvature scale becomes comparable to Ac, the photon dispersion relation 
is modified. The remarkable feature, however, is that this modification can induce 
a superluminal^ low-frequency phase velocity, i.e. the photon momentum becomes 
spacelike. 




Figure 1: Photons propagating in curved spacetime feel the curvature in the neighbourhood of 
their geodesic because they can become virtual e+e" pairs. 

At first, it appears that this must be incompatible with causality. However, 
as discussed in refs. [2-4], the relation of causahty with the "speed of hght" is far 
more subtle. For our purposes, we may provisionally consider causality to be the 
requirement that no signal may travel faster than the fundamental constant c defining 
local Lorentz invariancc. More precisely, we require that the wavefront velocity w^f, 
defined as the speed of propagation of a sharp-fronted wave pulse, should be less 
than, or equal to, c. Importantly, it may be shown [2,4,5] that v^i = t)ph(oo), the 
high-frequency limit of the phase velocity. In other words, causality is safe even if 
the low-frequency^ phase velocity i'ph(O) is superluminal provided the high-frequency 
limit does not exceed c. 

^In this paper, we use the term "superluminal" in the sense "greater than c". Apart from the 
occasional use of c in the text for clarity, wc set c = 1 throughout. Also, in our conventions, the 
metric of flat space is ry = diag (1,-1,— 1,-1) and the Riemann tensor is R^vaX = da^xv • 

^The term "low frequency" in this context requires some clarification. We work throughout in 
the WKB short wavelength approximation uj ^ R}/'^ and in the limit of weak curvature R <C 
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This appears to remove the potential paradox associated with a superluminal 
fph(O). However, a crucial constraint is imposed by the Kramers-Kronig dispersion 
relation^ (see, e.g. ref. [6], chpt. 10.8) for the refractive index, viz. 

Ren(oo) - Ren(O) = / — lmn{cu) . (1.1) 

Jo ^ 

where Iien{uj) = l/vpii{uj). The positivity of lmn{uj), which is true for an absorptive 
medium and is more generally a consequence of unitarity in QFT, then implies that 
Ren(cxD) < Ren(O), i.e. fph(c)o) > fph(O). So, given the validity of the KK dispersion 
relation, a superluminal Wph(O) would imply a superluminal wavefront velocity v^f = 
Vphioo) with the consequent violation of causality. 

We are therefore left with three main options [4], each of which would have 
dramatic consequences for our established ideas about quantum field theory: 



Option (1) The wavefront speed of light f^f > 1 and the physical light 
cones lie outside the geometric null cones of the curved spacetime, in 
apparent violation of causality. 



It should be noted, however, that while this would certainly violate causality for 
theories in Minkowski spacetime, it could still be possible for causality to be preserved 
in curved spacetime if the effective metric characterizing the physical light cones 
defined by i;„f nevertheless allow the existence of a global timelike Killing vector field. 
This possible loophole exploits the general relativity notion of "stable causality" [8,9] 
and is discussed further in ref. [2]. 



Option (2) Curved spacetime may behave as an optical medium ex- 
hibiting gain, i.e. Imn(co') < 0. 



This possibility was explored in the context of A-systems in atomic physics in ref. [4] , 
where laser-atom interactions can induce gain, giving rise to a negative Imn(cc;) and 
superluminal low-frequency phase velocities while preserving v^f = 1 and the KK 



where i? is a characteristic curvature of the background (which can also include derivatives of the 
curvature) and m is the electron mass. The frequency enters in the dimensionless ration cu'^R/m'^ 
and when we talk about "low" and "high" frequency we really mean small and large values of this 
dimensionless parameter. 

^Note that we are using "dispersion relation" in two different senses here. For clarity, we will 



always refer to eq.(Ll) explicitly as the Kramers-Kronig or KK dispersion relation to distinguish 
it from the use of the term dispersion relation to describe the frequency dependence of the photon 
light-cone. 
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dispersion relation. However, the problem in extending this idea to QFT is that 
the optical theorem, itself a consequence of unitarity, identifies the imaginary part 
of forward scattering amplitudes with the total cross section. Here, Im?7,(a;) should 
be proportional to the cross section for e'^e~ pair creation and therefore positive. A 
negative Imn(a;) would appear to violate unitarity. 

Option (3) The Kramers-Kronig dispersion relation ( p. . 1|) is itself vio- 
lated. Note, however, that this relation only relies on the analyticity of 
n{uj) in the upper-half plane, which is usually considered to be a direct 
consequence of an apparently fundamental axiom of local quantum field 
theory, viz. micro -causality. 

Micro-causality in QFT is the requirement that the expectation value of the com- 
mutator of field operators (0|[y4(x), y4(?/)]|0) vanishes when x and y are spacelike 
separated. While this appears to be a clear statement of what we would understand 
by causality at the quantum level, in fact its primary role in conventional QFT is 
as a necessary condition for Lorentz invariance of the S'-matrix (see e.g. ref. [6], 
chpts. 5.1, 3.5). Since QFT in curved spacetime is only locally, and not globally, 
Lorentz invariant, it is just possible there is a loophole here allowing violation of 
micro-causality in curved spacetime QFT. 

Despite these various caveats, unitarity, micro-causality, the identification of 
light cones with geometric null cones and causality itself are all such fundamental 
elements of local relativistic QFT that any one of these options would represent a 
major surprise and pose a severe challenge to established wisdom. Nonetheless, it 
appears that at least one has to be true. 

To understand how QED in curved spacetime is reconciled with causality, it is 
therefore necessary to perform an explicit calculation to determine the full frequency 
dependence of the refractive index n{uj) in curved spacetime. This is the technical 
problem which we solve in this paper. The remarkable result is that QED chooses 
option (3), viz. analyticity is violated in curved spacetime. We find that in the 
high-frequency limit, the phase velocity always approaches c, so we determine Wwf = 
1. Moreover, we are able to confirm that where the background gravitational field 
induces pair creation, 7 e~^e~, lmn{u) is indeed positive as required by unitarity. 
However, the refractive index n{uj) is not analytic in the upper half-plane, and the 
KK dispersion relation is modified accordingly. One might think that this implies a 
violation of microcausality, however, there is a caveat in this line of argument which 
requires a more ambitious off-shell calculation to settle definitively [7]. 
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In order to establish this result, we have had to apply radically new techniques to 
the analysis of the vacuum polarization for QED in curved spacetime. The original 
Drummond-Hathrell analysis was based on the low-energy, O^R/m?) effective action 
for QED in a curved background, 

C = -\f^uF^^ + [aRF^^F^^' + hR^,F^'^F\ + cR^^^^F^'' F^^'^ +■■■ . (1.2) 

derived using conventional heat-kernel or proper-time techniques (see, for example, 
[10-14]. A geometric optics, or eikonal, analysis applied to this action determines the 
low-frequency limit of the phase velocity. Depending on the spacetime, the photon 
trajectory and its polarization, Vph(O) may be superluminal [1,15,16]. In subsequent 
work, the expansion of the effective action to all orders in derivatives, but still at 
0{R/m'^), was evaluated and apphed to the photon dispersion relation [11,12,17, 
18]. However, as emphasized already in refs. [2,3,18], the derivative expansion is 
inadequate to find the high-frequency behaviour of the phase velocity. The reason is 
that the frequency uj appears in the on-shell vacuum polarization tensor only in the 
dimensionless ratio uj^R/m'^. The high-frequency limit depends non-pcrturbatively 
on this parameter'* and so is not accessible to an expansion truncated at first order 
in R/m^. 

In this paper, we instead use the world-line formalism which can be traced back to 
Feynman and Schwinger [19,20], and which has been extensively developed in recent 
years into a powerful tool for computing Green functions in QFT via path integrals 
for an appropriate 1-dim world-line sigma model. (For a review, see e.g. ref. [21].) 
The power of this technique in the present context is that it enables us to calculate the 
QED vacuum polarization non-perturbatively in the frequency parameter uj^Rjnv' 
using saddle-point techniques. Moreover, the world-line sigma model provides an 
extremely geometric interpretation of the calculation of the quantum corrections to 
the vacuum polarization. In particular, we are able to give a very direct interpretation 
of the origin of the Kramers- Kronig violating poles in niuS) in terms of the general 
relativistic theory of null congruences and the relation of geodesic focusing to the 
Weyl and Ricci curvatures via the Raychoudhuri equations. 

A further key insight is that to leading order in i?/m^, but still exact in cj^i?/m'*, 
the relevant tidal effects of the curvature on photon propagation are encoded in the 

^Notice that here we also include derivatives of the curvature in the generic symbol "i?". In 
fact, in ref. [18], the 0{R/rn?) contribution to the on-shell vacuum polarization was determined in 
the form n(a;) ~ ■:^f{^l ■ D)R, where £ ■ DR represents the variation of the curvature along the 
geodesic with tangent vector ^'^ and the function / is a form-factor determined from the effective 
action. This behaviour, where the vacuum polarization depends on the curvature through its 
variation 9„i?, where u is a light-cone coordinate adapted to the photon's original null geodesic, is 
reflected in the form of the Penrose limit for general curved spacetimes: see Section 7. 
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Penrose plane-wave limit [22, 23] of the spacetime expanded about the original null 
geodesic traced by the photon. This is a huge simplification, since it reduces the 
problem of studying photon propagation in an arbitrary background to the much 
more tractable case of a plane wave. In fact, the Penrose limit is ideally suited to 
this physical problem. As shown in ref. [24], where the relation with null Fermi 
normal coordinates is explained, it can be extended into a systematic expansion 
in a scaling parameter which for our problem is identified as R/rri^. The Penrose 
expansion therefore provides us with a systematic way to go beyond leading order in 
curvature. 

The paper is organized as follows. In Section 2, we introduce the world-line 
formalism and set up the geometric sigma model and eikonal approximation. The 
relation of the Penrose limit to the R/m^ expansion is then explained in detail, 
complemented by a power-counting analysis in the appendix. The geometry of null 
congruences is introduced in Section 3, together with the simplified symmetric plane 
wave background in which we perform our detailed calculation of the refractive index. 
This calculation, which is the heart of the paper, is presented in Section 4. The 
interpretation of the result for the refractive index is given in Section 5, where we 
plot the frequency dependence of n{(jj) and prove that asymptotically fph(c^) 1- 
We also explain exactly how the existence of conjugate points in a null congruence 
leads to zero modes in the sigma model partition function, which in turn produces 
the KK-violating poles in n{uj) in the upper half-plane. The implications for micro- 
causality are described in Section 6. Finally, in Section 7 we make some further 
remarks on the generality of our results for arbitrary background spacetimes before 
summarizing our conclusions in Section 8. 



2. The World-Line Formalism 




Figure 2: The loop x'^{t) with insertions of photon vertex operators at t\ and T2- 
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In the world-line formalism for scalar QED^ the 1-loop vacuum polarization is 
given by 



a 



I 

Jo 



oo 



dT 



■T 



n 



1-loop 



dndT2Z{v:jxin)]v^,,Mr2)]) . 



(2.1) 



The loop with the photon insertions is illustrated in Figure (|]). The expectation 
value is calculated in the one-dimensional world-line sigma model involving periodic 
fields x^(r) = x^(r + T) with an action 



where m is the mass of the (scalar) electron and we work in Minkowski signature 
in both spacetime and on the world-line.^ The factor Z is the partition function of 
the world-line sigma model relative to flat space.'' It is an important detail of our 
calculation that Z will depend implicitly on u and the insertion points ri and T2. 

The vertex operators have the form 



where A^{x) is the gauge connection of a photon propagating with momentum k 
and polarization vector e. At the one-loop level, we can impose the tree-level on- 
shell conditions for the gauge field. This means D^F^^ = along with the gauge 
condition D^A^ = 0. In curved spacetime, the photon gauge field is not exactly that 
of a plane wave due to the effects of curvature and in general it would be impossible 
to solve for the on-shell vertex operator. However, we will work in the WKB, or 
short wavelength, approximation which is valid when to ^ i?^/^.^ This is the limit 

^Since all the conceptual issues we address are the same for scalars and spinors, for simplicity 
we perform explicit calculations for scalar QED in this paper. The generalization of the world-line 
formalism to spinor QED is straightforward and involves the addition of a further, Grassmann, field 
in the path integral. For ease of language, we still use the terms electron and positron to describe 
the scalar particles. 

^This will require some appropriate ie prescription. In particular, the T integration contour 
should lie just below the real axis to ensure that the integral converges at infinity. 

''in general, one has to introduce ghost fields to take account of the non-trivial measure for the 
fields. 



in curved spacetime [25-29]. However, in our calculation where we work to leading order in R/m? 
in a special set of coordinates the determinant factor is 1 to leading order. 

^It is important to understand that this notion of high frequency still allows one to expand the 
effective action in powers of lo because this latter is actually a function of the dimensionless ratio 
bj'^R/m'^ which can be small. 




(2.2) 
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of geometric optics where Af^{x) is approximated by a rapidly varying exponential 
times a much more slowly varying polarization. Systematically, we have 

A^{x) = (e^ix) + co-'B^ix) + ■■■ )e^-®(^) . (2.4) 

We will need the expressions for the leading order pieces and e. This will necessitate 
solving the on-shell conditions to the first two non-trivial orders in the expansion in 
R^/'^ /uj. To leading order, the wave-vector = coi^, where = 9^9 is a null vector 
(or more properly a null 1-form) satisfying the eikonal equation, 

£■£ = g^^'df.ed^Q = . (2.5) 

A solution of the eikonal equation determines a family or congruence of null geodesies 
in the following way.^ The contravariant vector field 

r(x) = d^'Qix) , (2.6) 

is the tangent vector to the null geodesic in the congruence passing through the point 
x'^. In the particle interpretation, = ui^ is the momentum of a photon travelling 
along the geodesic through that particular point. It will turn out that the behaviour 
of the congruence will have a crucial role to play in the resulting behaviour of the 
refractive index. The general relativistic theory of null congruences is considered in 
detail in Section 3. 



Now we turn to the polarization vector. To leading order in the WKB approxima- 
tion, this is simply orthogonal to i, i.e. e-i = 0. Notice that this does not determine 
the overall normalization of e, the scalar amplitude, which will be a space-dependent 
function in general. It is useful to split e'^ = Ai'^, where is unit normalized. At the 
next order, the WKB approximation requires that is parallel transported along 
the geodesies: 

i-De'' = 0. (2.7) 
The remaining part, the scalar amplitude A, satisfies 

i- DlogA = --D -i . (2.8) 

Eqs. ( |2.7|) and ( p.8|) are equivalent to 

i- De^ = -^e^'D -i . (2.9) 

Since the polarization vector is defined up to an additive amount of k, there are two 
linearly independent polarizations ei{x), i = 1,2. 

^Thc congruence is not, in general, unique due to existence of integration constants. Later we 
will find that our results are independent of these integration constants. 
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Since there are two polarization states, the one-loop vacuum polarization is ac- 
tually a 2 X 2 matrix 



In order for this to be properly defined we must specify how to deal with the zero 
mode of x^{t) in the world-line sigma model. Two distinct - but ultimately equiv- 
alent - methods for dealing with the zero mode have been proposed in the litera- 
ture [25-29]. In the first, the position of one particular point on the loop is defined 
as the zero mode, while in the other, the "string inspired" definition, the zero mode 
is defined as the average position of the loop: 



We will use this latter definition since it leads to a much simpler formalism. Since 
we are effectively calculating an on-shell term in an effective action, the integral over 
the zero mode is simply excluded from the functional integral. In other words, our 
world-line sigma model does not include an integral over Xq which one should think of 
as being a fixed point in spacetime. Since in curved spacetime there is in general no 
translational symmetry, the one-loop correction n^^''°°^(a;o) will depend explicitly on 
Xq. We will always choose coordinates for which Xq = 0, in which case we implicitly 
impose the constraint 



on the sigma model fields. The advantage of using the string inspired method is that 
there is translational symmetry on the world-line loop. This allows us to fix ri = 0. 
We will then take ti = ^T, < ^ < 1, and replace the two integrals over ti and T2 
by a single integral over the variable ^. 

A key ingredient in our analysis is that in the limit of weak curvature R <^ 
m^, the sigma model based on the general metric Qf^i, can be approximated by the 
metric in a cylindrical neighbourhood of the geodesic in the null congruence that 
passes through Xq = 0. We will call this particular geodesic 7. The metric in 
the neighbourhood of 7 arises in a very particular way known as the Penrose limit 
[22]. Exactly how this limit arises is rather remarkable and means that the vacuum 
polarization and refractive index is only sensitive to the Penrose limit of the original 
metric. It should be noted that the Penrose hmit captures the global behaviour of 
the original metric all the way along the geodesic 7. 




(2.10) 




(2.11) 




(2.12) 
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Now notice that the exponential pieces of the vertex operators in ( |2.1|) act as 
source terms and so the complete action including these is 



S = -T + — J dTg^,{x)x^x''-uQ[x{O]+uje[x{0)]. (2.13) 

Here, we have scaled r — >■ Tr and then T T/m^, so that r runs from to 1. T is 
now dimensionless and plays the role of a conventional coupling constant. In 

fact, the effective coupling constant is actually the dimensionless ratio R/m?, where 
i? is a typical curvature scale. So when R/rn^ is small we can perform a perturbative 
expansion in the world-line sigma model. As is usual in a perturbative analysis, it is 
useful to re-scale the "fields" x^{t) appropriately in order to remove the overall factor 
of /T. The coupling then re- appears in vertices. However, this re-scaling must be 
done in a clever way. The reason is that the classical saddle-point solution following 
from ( p.l3|) is not simply the constant configuration x^(r) = Xg = because the 
sources inject world-line momentum into, and out of, the system. It is not difficult 
to guess what the classical saddle-point solution will be because the classical equation 
of motion that follows from ( p.l3| ) is just the geodesic equation for x^ij) with delta- 
function sources at r = and ^. The solution consists of an electron and positron 
pair produced at r = = 1 which propagate along the photon geodesic 7, with the 
electron going from r = to r = ,^ and the positron from r = 1 to r = ^, before 
annihilating at r = ,^ back into a photon which then continues along 7. In other 
words, the classical loop is squashed onto the geodesic 7 as illustrated in Figure @). 
We will find the explicit solution for this classical loop shortly. 











e- 





Figure 3: The classical saddle point solution consists of a squashed loop which follows the geodesic 
7. The length of the loop is ~ Lo/m? which represents a potentially interesting UV-IR mixing effect. 



As we have said, the fact that there is a non-trivial classical solution around 
which the perturbative expansion is performed means that the re-scaling of the fields 
must be done in an appropriate way. The problem is solved by choosing from the 
outset a set of coordinates which are adapted to the null congruence containing 7. 
These coordinates (u, 6, y"), a = 1, 2, are known as Rosen coordinates . They include 
two null coordinates: m, the affine parameter along the geodesies and 6, the solution 
of the eikonal equation ( ^.51 ). As explained in ref. [23], the full metric g^^ around 7 
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can always be brought into the form 

ds^ = 2du dQ - C{u, e, Y'')dQ^ - 2Ca{u, 0, Y'')dY'' dO - Cab{u, 6, Y'')dY^ dY^ . 

(2.14) 

It is manifest that dQ is a null 1-form. The null congruence has a simple description 
as the curves (m, 9o, Yq) for fixed values of the transverse coordinates (Go, Yq). The 
geodesic 7 is the particular member (m, 0, 0, 0). It should not be surprising that the 
Rosen coordinates are singular at the caustics of the congruence. These are points 
where members of the congruence intersect and will be described in detail in the next 
section. 



With the form ( p.l4| ) of the metric, one finds that the classical equations of 



motion of the sigma model action ( |2.13| ) have a solution with = O = where 
u{t) satisfies 

^ = -^^(r-0 + ^5(r). (2.15) 

More general solutions with constant but non-vanishing (B, Y"") are ruled out by the 
constraint ( p.l2|) . The solution of ( p.l5|) is 



~( \ ^/2a;T(l-Or/m2 0<r<e 

uir) = —uq + < (2.16) 
[2ujT^{l-T)/m^ e<r< 1 . 



where the constant 

Uo = ujT^{l-^)/m^ (2.17) 



ensures that the constraint ( p.l2|) is satisfied. The solution describes a loop which 



is squashed down onto the geodesic 7 as illustrated in Figure (|^). The electron and 
positron have to move with different world-line velocities in order to accommodate 
the fact that in general ^ is not equal to |. In Section 5, we explain how for particular 
values of T there are more general classical saddle-point solutions which are consistent 
with ( |2.12| ). However, the solution we have described is the only one that exists for 



generic values of T. 

What is intriguing about this picture is that the classical loop, which has an affine 
parameter length proportional to L ~ oj/rn?, actually gets bigger as the frequency 
is increased. The reason is that higher frequency leads to bigger impulses and hence 
longer loops. This is an interesting example of the kind of UV-IR mixing that is 
seen in other contexts, such as non-commutative field theories or high energy string 
scattering. Whether the occurrence here is hinting at something deeper deserves 
to be investigated in more detail. However, what it will mean is that the higher 
frequencies will probe global aspects of the spacetime rather than shorter distance 
scales as our intuition might have suggested. 
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Now that we have defined the Rosen coordinates and found the classical saddle- 
point solution, we are in a position to set up the perturbative expansion. The idea is 
to scale the transverse coordinates B and in order to remove the factor of m} jT 
in front of the action. The affine coordinate m, on the other hand, will be left alone 
since the classical solution u[t^ is by definition of zeroth order in perturbation theory. 
The appropriate scalings are precisely those needed to define the Penrose limit [22] - 
in particular we closely follow the discussion in [23]. The Penrose limit involves first 
a boost 

(u, 6, Y"") — > (A-^M, AG, y") , (2.18) 
where A = T^l"^ jm^ and then a uniform re-scaling of the coordinates 



(AM,Ae,Ar" 



(2.19) 



As argued above, it is important that the null coordinate along the geodesic u is not 
affected by the combination of the boost and re-scaling; indeed, overall 



After these re-scalings, the sigma model action ( p.l3|) becomes 



(2.20) 



S 



T + 



4 



dr 



2uQ-\^c{u, A2e,Ay'')e^ 



2Aa(M, A^e, AF^)F'* - CaAu. A^O, Ar")F" Y 



(2.21) 



^e(o + ^e(o). 



In the limit R -C m^, we expand in powers of A = T^l"^ jm and ignore terms of 0{X): 



S 



dr 



uT 



ujT 



2ue-Cab{u,o,o)Y''Y'' - e(o + ^e(o) + -- - . (2.22) 



The leading order piece is precisely the Penrose limit of the original metric in Rosen 
coordinates. Notice that we must keep the source terms because the combination 
ujT/m?, or more precisely the dimensionless ratio toR^^'^ /m^, can be large. However, 
there is a further simplifying feature: once we have shifted the "field" about the clas- 
sical solution u{t) uIt) +u{t), it is clear that there are no Feynman graphs with- 
out external B lines that involve the vertices d^Cab{'u,0,0)u"'Y°'Y^, n > 1; hence, 
we can simply replace Cab{u + u, 0, 0) consistently with the background expression 
Cab{u, 0, 0). This means that the resulting sigma model is Gaussian to leading order 
in R/m?: 



dr 



2Me-Cafe(u,o,o)r"r' 



dTCab{u,^,m"y' , (2.23) 
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where finally we have dropped the iiQ piece since it is just the same as in fiat space 
and the functional integral is normalized relative to fiat space. This means that 
all the non-trivial curvature dependence lies in the Y°- subspace transverse to the 
geodesic. 

It turns out that the Rosen coordinates are actually not the most convenient co- 
ordinates with which to perform explicit calculations. For this, we prefer Brinkmann 
coordinates {u, v, y*). To define these, we first introduce a "zweibein" in the subspace 
of the Y"": 

CM = 5,,E\{u)E\{u) , (2.24) 
with inverse E°'i. This quantity is subject to the condition that 

dE 

^ = (2.25) 

^ du ^ ^ ^ 

is a symmetric matrix. Then the affine coordinate u is common to both systems, 
while 

yi ^ ^ijra ^ v = Q + -'^EW^'Y^ . (2.26) 

2 du 

Notice that the Brinkmann coordinates are homogeneous under the scaling ( p. 20 ): 

{u,v,y')-^{uA^v,\y') . (2.27) 



In Brinkmann coordinates, the metric takes the form 

ds^ = 2du dv + h,j{u)y' du^ - dy^ , (2.28) 
where the quadratic form is 

K,{u) = -^E^j ■ (2.29) 

We have introduced these coordinates at the level of the Penrose limit. However, 
they have a more general definition for an arbitrary metric and geodesic. They are 
in fact Fermi normal coordinates. These are "normal" in the same sense as the 
more common Riemann normal coordinates, but in this case they are associated to 
the geodesic curve 7 rather than to a single point. This description of Brinkmann 
coordinates as Fermi normal coordinates and their relation to Rosen coordinates and 
the Penrose limit is described in detail in ref. [24]. In particular, this reference gives 

"'^"An alternative proof of this result which relies only on conventional power counting arguments 
and does not rely on any a priori knowledge of the Penrose limit is provided in Appendix A. 
"'^^Notice that i and j are raised and lowered in this Euclidean 2d subspace by 6ij and not — 
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the A expansion of the metric in null Fermi normal coordinates to (9(A^). To 0{X) 
this is 



ds'^ = 2du dv — RiujuU^ V'' dv^ — dy 



i2 



+ A 



- 2Ruiuvy'v du^ - ^RuijkVY du dy^ - \Ruiur,ky'y^y^ du^ + 0{\^) , 



(2.30) 



which is consistent with ( p.28| ) since Riuju = —hij for a plane wave. It is worth 
pointing out that Brinkmann coordinates, unlike Rosen coordinates, are not singular 
at the caustics of the null congruence. One can say that Fermi normal coordinates 
(Brinkmann coordinates) are naturally associated to a single geodesic 7 whereas 
Rosen coordinates are naturally associated to a congruence containing 7. 



In Brinkmann coordinates, the Gaussian action ( p.23| ) for the transverse coordi- 
nates becomes 



^^iiVV'' 

r=5 2m? ^ 



T = 



(2.31) 



where the world-line velocity along the loop is 

u(t) -- 



2luT{1- ^)/m^ 0<r<^ 



(2.32) 



Although ( |2.31| ) looks more complicated than ( p.23| ), it is actually more useful for 
explicit calculations. 



3. The Symmetric Plane Wave and Null Congruences 



The analysis above shows that photon propagation in a completely general curved 
spacetime is governed to one-loop order by the Penrose limit for the metric in a neigh- 
bourhood of the original null geodesic. The complete one-loop vacuum polarization 
and photon dispersion relation can therefore be determined without loss-of-generality 
by working in a plane wave background. 

In Section 7, we briefly discuss the Penrose limits of spacetimes of special physical 
interest, such as de Sitter and Schwarzschild, and see how known results for low- 
frequency photon propagation in these spacetimes are recovered as special properties 
of the Penrose limit. For the rest of this paper, however, we specialize to the simplest 
example of a plane wave - the symmetric plane wave [23]. In this background, we 
can evaluate the non-perturbative frequency dependence of the vacuum polarization 
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explicitly. In doing so, we discover many surprising features of the dispersion relation 
that will hold in general. 



The symmetric plane wave metric is given in Brinkmann coordinates by ( |2.28| ), 
with the restriction that hij is independent of u. This metric is locally symmetric 
in the sense that the Riemann tensor is covariantly constant, DxRfj^^pcr = 0, and can 
be realized as a homogeneous space G/H with isometry group G}^ With no loss 
of generality, we can choose a basis for the transverse coordinates in which hij is 
diagonal: 

h.jyy = al{y'f + al{y'f. (3.1) 

The sign of these coefficients plays a crucial role, so we allow the cTj themselves to be 
purely real or purely imaginary. 

For a general plane-wave metric, the only non-vanishing components of the Rie- 
mann tensor (up to symmetries) are 

Ruiuj — hij{u) . (3-2) 
So for the symmetric plane wave, we have simply 

Ruu = C^l + 0"2 5 o\ 

R - -^2 ^^-^^ 



and for the Weyl tensor. 



i=i 

The null energy condition, viz. T^^k'^k'^ > with k^ a null vector, reduces here to 
Tuu > 0, so from Einstein's equation we require Ruu = o"^ + (t| > 0. It follows that 
at least one of the cTj must be real (we will always choose this to be (Xi). Special 
choices of the CTj allow the symmetric plane wave to be either Ricci flat (cxi = ±icr2) 
or conformally flat (ai = ±cr2). The Ricci fiat case is the vacuum gravitational wave. 



While, as we saw in the last section, the original null geodesic 7 (with i = du) 
defines the classical solution in the world-line path integral, in order to evaluate the 
fiuctuations we also need the eikonal phase and wave-vector for deviations from 7 
itself. We therefore need to study the congruence of null geodesies in the neighbour- 
hood of 7 in Brinkmann coordinates. We first do this explicitly for the symmetric 

^^Notice that, contrary to the imphcation in ref. [4, 18], the condition that the Riemann tensor 
is covariantly constant only implies that the spacetime is locally symmetric, and not necessarily 
maximally symmetric [13,23]. A maximally symmetric space has R^,ypa — ]^^(3/^p5i/ct ~ 9ij,agvp) 
and does not have the required anisotropy for the vacuum polarization to modify the speed of light. 
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plane wave background, then explain how the key features are described in the gen- 
eral theory of null congruences. 



The geodesic equations for the symmetric plane wave (|2.28|) , (|3.1| ) are: 

u = , 

2 

v + 2uY,^hY = , (3.5) 

i=l 

We can therefore take u itself to be the affine parameter and, with the appropriate 
choice of boundary conditions, define the null congruence in the neighbourhood of, 
and including, 7 as: 

V = Q > CTj tanfcTjM + aAy^'^ , 

2 (3.6) 

?/* = cos(crjU + ttj) . 

The constants B and are nothing other than the Rosen coordinates for the sym- 
metric plane wave. In fact, in Rosen coordinates the symmetric plane wave metric 
is 

2 

ds"^ = 2du dQ-Y^ cos^{(JiU + ai)dY''^ . (3.7) 

i=l 

The integration constants can be thought of as redundancies in the definition of the 
null congruence and the associated Rosen coordinates; in particular, they determine 
the position of the caustics. Given this, we have 



E^'i = 5ia sec{aiU + at) , (3.8) 
flij = —SijCi tan(crjU + a^) 



and it is immediate that the eikonal phase is 

2 

2 



9(x) = f H — o"j tan(o"iM + a^)?/*^ . (3.9) 

■t=i 

The tangent vector to the congruence, defined as i'^ = g^'^d^Q, is therefore 
^ 2 

^ = du + -Y^ {fx- (tan2((Tin + a^) - l)y'^ d., - aitixn{aiU + ai)y'd,^ . (3.10) 



1=1 
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The polarization vectors are orthogonal to this tangent vector, I ■ Si = 0, and are 
further constrained by (|2.9| ). Solving (|2.71) for the normalized polarization (one-form) 
yields^^ 

Si = dy"^ + (Tj tan(crj'u + ai)y^du . (3-11) 

The scalar amplitude A is determined by the parallel transport equation ( p.8|) , from 
which we readily find (normalizing so that ^(0) = 1) 



^ = 11^/ . (3-12) 

V COS((TjM + ai) 
1=1 ' 



The null congruence in the symmetric plane wave background displays a number 
of features which play a crucial role in the analysis of the refractive index. They 
are best exhibited by considering the Raychoudhuri equation, which expresses the 
behaviour of the congruence in terms of the optical scalars, viz. the expansion 9, 
shear a and twist u). These are defined in terms of the covariant derivative of the 
tangent vector as [30]: 

= ^D/^ , 



UJ 



\D^^^DH--e^ , (3.13) 



The Raychoudhuri equations describe the variation of the optical scalars along the 
congruence: 

3.14 

(9„<T = -2^a- l^ol • 

(We will not need the equation for the twist.) Here, we have introduced the Newman- 
Penrose notation (see, e.g. ref. [30]) for the components of the Ricci and Weyl tensors: 
$00 = ^R^Jt^r, ^0 = Cf.p^Jf'rmPm''.^^ As demonstrated in refs. [31], the effect 
of vacuum polarization on low-frequency photon propagation is also governed by 
the two curvature scalars $oo and \l/o- Indeed, many interesting results such as 
the polarization sum rule and horizon theorem [31,32] are due directly to special 
properties of $00 and \l'o. As we now show, they also play a key role in the world- line 
formalism in determining the nature of the full dispersion relation. 



^^The one- form is exactly what appears in the vertex operator via ei^i^. 

^■^For the symmetric plane wave, the Newman-Penrose null tetrad basis n^, m^, comprises 
£ as in eq. (3.1C), n — dy, and m = "^(si +*£2)- The basis vectors satisfy £ ■ n — 1, m ■ fh — —1 and 



the metric can be expressed as g^j/ — 2{£(^^ny-^ — m(^^m^^Y Since eqs. (3^) and (3^) are the only 



non- vanishing components of the Ricci and Weyl tensors, it follows that $00 = ^Ruu = +^2 

^0 — ^(C„i„i — Cu2u2) — \{'^2 
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By its definition as a gradient field, it is clear that -D^^f^,] = so the null con- 
gruence is twist-free uj = 0. The remaining Raychoudhuri equations can then be 
rewritten as 

dui9 + a) = -{§ + af - $00 - |^o| , ^ 

(o.loj 

- a) = -(e-^)2- $00 + 1^01 . 

The effect of expansion and shear is easily visualized by the effect on a circular cross- 
section of the null congruence as the affine parameter u is varied: the expansion 6 
gives a uniform expansion whereas the shear a produces a squashing with expansion 
along one transverse axis and compression along the other. The combinations ^ ± o" 
therefore describe the focusing or defocusing of the null rays in the two orthogonal 
transverse axes. 

We can therefore divide the symmetric plane wave spacetimes into two classes, 
depending on the signs of $oo ± |^o|- A Type I spacetime, where $00 ± |^o| are 
both positive, has focusing in both directions, whereas Type II, where $00 =t \E'o 
have opposite signs, has one focusing and one defocusing direction. Note, however, 
that there is no "Type III" with both directions defocusing, since the null-energy 
condition requires $00 > 0. 

For the symmetric plane wave, the focusing or defocusing of the geodesies is 
controlled by eq.( p.6| ), = cos(criM + Oj). Type I therefore corresponds to cxi and 
(72 both real, whereas in Type II, ai is real and (T2 is pure imaginary. The behaviour 
of the congruence in these two cases is illustrated in Figure (^). 





Figure 4: (a) Type I null congruence with the special choice (Ti ~ 0-2 and ai = 02 so that 
the caustics in both directions coincide as focal points, (b) Type II null congruence showing one 
focusing and one defocusing direction. 



To see this explicitly in terms of the Raychoudhuri equations, note first that the 
curvature scalars $00 ~ l^ol = <^i, ^00 + l^ol = c^i are simply the eigenvalues of hij. 
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The optical scalars are 

d = — (cTi tan(criM + ai) + a2 tan(cr2U + 02)) , 

1 2 (3.16) 
0" = - (cTi tan((Titi + Oi) — (T2 tan(cr2M + 02)) 

and we easily verify 

1 ^ (3.17) 
du<y = -2ea+-{al-al) . 

It is clear that provided the geodesies are complete, those in a focusing direction 
will eventually cross. In the symmetric plane wave example, with ?/* = cos(crjM + 
ttj), these "caustics" occur when the affine parameter cTjU = 7r(n + |) — aj, n G Z. At 
a caustic, the amplitude factor A in (|3.12|) diverges and correspondingly the Rosen 
coordinates are not well defined. In fact, the existence of conjugate points, i.e. points 
p and g on a geodesic 7 that can be joined by geodesies infinitesimally close to 7, is 
generic in spacetimes satisfying the null energy condition. The result is summarized 
in the following theorem [8,33]: 

Theorem: if a spacetime satisfies the "null generic condition" [i.e. every 
null geodesic has at least one point where either 

R^jn'' ^ or lixC,^pu[Jr]^'t ^ 0, (3.18) 

or equivalently $00 7^ or \E'o 7^ 0) and the null energy condition, then 
every complete null geodesic possesses a pair of conjugate points. 

The existence of conjugate points will turn out to be crucial in the world-line 
sigma model formalism. It means that for certain values of T (for a given a;), such 
that u = ±uo are conjugate points, in the Penrose limit around the geodesic, there 
exists a family of classical solutions corresponding to the different geodesic paths 
between the conjugate points. This implies the existence of zero modes which, as 
explained in Section 5, ultimately controls the location of singularities of the refrac- 
tive index in the complex u plane and is the key to understanding the violation of 
the conventional Kramers-Kronig dispersion relation and the fate of micro-causality. 

"'^^This does not necessarily mean that the conjugate points are joined by more than one actual 
geodesic, only that an infinitesimal deformation of 7 exists. Later we shall see that the existence 
of conjugate points relies on the existence of zero modes of a linear problem. Conversely, the 
existence of a geodesic other than 7 joining p and q does not necessarily mean that p and q are 
conjugate [8,33]. 

^^Whether these deformed geodesies become actual geodesies is the question as to whether they 
lift from the Penrose limit to the full metric. 
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4. World-line Calculation of the Refractive Index 



In this section, we calculate the vacuum polarization and refractive index explicitly 
for a symmetric plane wave. As we mentioned at the end of Section 2, the explicit 
calculations are best performed in Brinkmann coordinates. We will need the expres- 
sions for G and Ei for the symmetric plane wave background: these are in eqs. (|3.9|) , 
( |3.11|) and ( p.l2| ). From these, we have the following explicit expression for the vertex 



operator^'' 



cos a,- 



X expzcu 



1 ^ 

f + - CTj tan(crjM + aj)y^'^ 



(4.1) 



The Gaussian action for the transverse coordinates, including the source terms (|2.31| ), 
is 

i=i ^ (4.2) 
- ( tan(o-iUo + aj)?/'(0^ + tan(criMo - ai)yX0Y) | . 

Notice that the fluctuations are completely decoupled. The measure for the field 
x^{t) is covariant and so includes the factor a/— det (7[a;(r)] which can be exponenti- 
ated by introducing appropriate ghosts [25-29]. However, in Brinkmann coordinates 
after the re-scaling ( ^.27] ), det g = —1 + 0{X) and so to leading order in R/rn? the 
determinant factor is simply 1 and so plays no role. The same conclusion would not 
be true in Rosen coordinates. 



The fluctuations satisfy the eigenvalue equation 

y' + u^a}y' -^{i&n{aiUo + ai)5{T -C) + ^a.n{aiUo-ai)5{T))y' = Xy'-C , (4.3) 

where C is the Lagrange multiplier that is determined by imposing the constraint 
£dTy' = 0. Now we see the utility of the Brinkmann coordinates, because the 
equation (|4.3| ) is just that of a simple harmonic oscillator and the non-trivial aspects 
of the problem lie solely in the matching conditions at r = and r = ^. On the 
contrary, in Rosen coordinates the eigenvalue equation has hypergeometric solutions 

-'^ "^Notice that at leading order in R/m? we are at liberty to replace u(r) by its classical value 
u{t). The argument is identical to the one given in Section 2. 
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and is not so straightforward to deal with. Consequently, we search for a solution in 
the form 

... (AiCos{uJir) + Bism{uJiT)-C/ujj 0<t<^ 
y (t) ^ \ (4-4) 
\^A2 cos{u!2t) + B2 sin(a;2T) — C/a;| ^ < t < 1 . 

where ojf = Au;^T^af{l - ifjm^ - A and a;| = AiJ^T^ j - X. The matching 
conditions at r = 0(= 1) and t — ^ are the continuity of and the jumps 

Ay*(0) = ^^tan(a,«o - ai)y%Q) , 
= 5-^tan((7iMo + ai)y\C) ■ 

These conditions, along with dTy^(T) = 0, determine the five unknowns Ai, Bi 
and C. A solution is only possible if A satisfies a characteristic equation ^(A) = 0. 
When T{X) is suitably normalized, the determinant of the fluctuation operator is 
given by J^(0). This leads to the remarkably simple formula for the determinant 
factor relative to flat space: 



■^j^ y cos^ tti sm pj cos Pi 

where 

ft = ^Mfli . (4.7) 

Notice that Z — > 1 in the flat space limit CTj — > 0. 

The remaining correlation function piece is determined by the Green function 

G,j{ry)^{yW{T')) . (4.8) 

It is clear that this is diagonal in the polarization indices, where the diagonal com- 
ponents are the solution of the equation 

dl + u{r fai 5-^ ( tan(criMo + ai)6{r - ^) 

(4.9) 

+ tan((7iiio - ai)5{T)) Guir, r') = -2i5{T - t') - C , 

We can find this by a brute force solution similar to that above, imposing boundary 
conditions so that Gn(r, r') is continuous at r = 0, r' and ^ and its derivative jumps 
by the appropriate amounts at r = 0, r' and ^. As before, C is determined by 
imposing dr Gii{T,T') — 0. The solutions themselves are not very illuminating 
and so we do not write them down here. The Green function leads to the following 
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remarkably simple formula for^^ 



^ ^1 V cos(c7iM(r) + ai) cos(cr/M(r') + ai) 
X ^drdr' + aitaxi{aiu{T) + ai)u{T)dr' + aita,n{aiu{T') + ai)u{T')dr 
+ af iaji{aiu{T) + a^) tan(cri'u(r') + ai)-u(r)'u(r') Gii(r, r') 



t=^,t'=0 



T I sin Pi cos A >' fj- V cos(A + cos(A - ai) " 

(4.10) 

The i here is crucial and appears because we are working in Minkowski signature. 

Putting all these pieces together, the final result for the one-loop correction to 
the vacuum polarization is 

"OO jrp pi 



am'' p dT . f\, \, A 1^ / Pf 



2n Jo T' Jo [ sinAcosA V sin^ AcosA 

(4.11) 

The term 5(^)Z(A)/(cosAcos A) has been replaced by 1 since 



2 



\im Z{Pi)Y\ J ^ = 1 . (4.12) 

^^'^ 11 y cos(A + a;) cos(A - az) ^ ' 

It is remarkable that the result for the vacuum polarization is independent of Oj 
so the ambiguity in the choice of the null congruence has no effect on the final 
result. It is especially noteworthy that the divergences of the vertex operators due 
the singularities of the scalar amplitude at the caustics of the null congruence are 
completely removed by quantum effects. 

The mass-shell conditions for the two polarization states are modified by the 
one-loop correction to 

l(a;^-^^) + nl:'-» = 0. (4.13) 



^*The limits r — > ^ and r' — > have to be taken after the derivatives have been evaluated due to 
implicit dependence on ^. 
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The phase velocities are Vph = ^/\k\ then and hence the refractive indices for the 
two velocity eigenstates are 



1^1 + 2nj:'°°PH 1 ^, 



^ + 27ru;2 ^2 j 1 n y sln^ A cos A 



(4.14) 



to order a. In particular, notice that the polarization vectors Si correspond directly 
to the two velocity eigenstates. 



5. Analysis and Interpretation 



The first remark is that the expression for the refractive indices (|4.14|) is completely 
UV safe since the term in curly brackets behaves as for small T. This can be 
traced to the fact that we have imposed the tree-level on-shell condition on the 
photon momentum. 

As we proceed, it is useful to have in mind the behaviour of the refractive index 
in a simple model of a dissipative dielectric medium with a single absorption band.^^ 
This is modelled by an electric permittivity of the form 

ji 

e{uj) = 1 - . (5.1) 

where uo is the resonant frequency and 7 is the width. For weak coupling, 

n{uj) = ^) = l- ^ ""'Z^. +■■■ . (5.2) 



Written in the same form as ( [4.14[ ) as a T integral, we have 



n{uj) = l ^ dTe-'^e-^^'^/^^^o)sm{iuT/iUo) . (5.3) 

^i^o Jo 

In this case, the T integral is perfectly well defined without the need for an ie 
prescription. The real and imaginary parts of n{uj) — 1 are sketched in Figure (|^). 
At low frequencies the phase velocity is subluminal. At frequencies ~ cjq the 
imaginary part of n{uj) has an absorption peak and the phase velocity changes over 
to being superluminal. At high frequencies, the phase velocity approaches 1 as 1/a;^. 
It is important to emphasize that the superluminal phase velocity at high frequencies 
is not associated with a violation of causality since asymptotically it approaches c. 

-^^This simple model forms the basis of many textbook discussions; for example, see Jackson [34], 
chpt. 7.10. 
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5.1 The analytic structure of the integrand 



When we compare our result ( [4.14D to the simple model of a dissipative medium 



]3|), the most striking difference is the existence of singularities in the integrand. 
When (Tj is real, the integrand ( [4.14[ ) has branch point singularities on the positive 
real axis at 



TTm?n 



T = , n = l,2,... (5.4) 

and the T integral must be properly defined in order to have a finite result. The 
correct procedure is to take the contour to lie just below the real axis. It is significant 
that these singularities arise from zeros of the fiuctuation determinant ( |4.6| ) and have 
a natural interpretation in terms of zero modes, viz. non-trivial solutions of ( |4.^j| ) 
with zero eigenvalue A = 0. For the special case when i = \ these zero modes are 
particularly simple: u = u{t) as in ( p.l6|) and 

y'ij) = sin(2n7rr) . (5.5) 



The expression for t>(r) is then completely determined by solving the geodesic equa- 
tion 

V + J2 {^^^IvY + ^y^' sec2(a,« + a,) (5(r - - 5(r)) } = . (5.6) 

i=\ 

These solutions are therefore associated with geodesies that are arbitrarily close to 
7 that intersect 7 at both u = ±mo and for > 1 at points in between. In other 
words, u = ±uo are conjugate points on the geodesic 7. The n = 1 and n = 2 zero 
modes are illustrated in Figure (^. 
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Figure 6: The n = 1 (red) and n — 2 (blue) zero modes for ^ = i. The points u = ±uo are 
conjugate points for 7. 



For generic ^ the solutions are more complicated: again u = u{t) as in ( p.l6| ) 
while 

. ( Ai sm{aiu{T)) + Bi cos((JiM(r)) < r < ^ 

y \t) = < ^ ^ (5.7) 

I A2 sin(aiu(r)) + i?2 cos((Tin(r)) ^ < r < 1 



and f (r) solves ( |5.6| ). Imposing the continuity of ?/* and the conditions ( [4 .51 ) implies 
that for n odd we must have 

Ai=A2 = ^^^5itana, , B2 = "J^^^i ' (5-8) 

and so there is only a single zero mode. For n even, however, there are two zero 
modes since there are only two conditions on the four constants: 

B^ = B2, (l-0A + e^2 = 5itanai . (5.9) 

Once again these solutions are formed from portions of two inequivalent geodesies 
which intersect at x^{0) and x'^{^): in other words, the points x'^(O) and x^{^) are 
conjugate, however, for Oj 7^ the conjugate points do not generally lie on 7. An 
example of the first zero mode is illustrated in Figure (|^. 

Similarly the singularities on the imaginary T axis that arise when one of the a", 
is imaginary can be understood in terms of these zero modes, but with imaginary 
affine parameter u iu. It is tempting to think that these solutions in imaginary 
affine parameter can be associated with world-line instanton solutions. These kinds 
of instanton in the world-line sigma model (not to be confused with instantons in the 
original field theory) have been discussed in the literature in the context of non-trivial 
electromagnetic backgrounds and can be used to describe Schwinger pair creation in 
that context [35,36]. Here, we shall shortly see why this interpretation is exactly 
right, since the singularities on the imaginary axis determine the imaginary part of 
the refractive index which describes the dissipative nature of the propagation that 
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Figure 7: The n = 1 zero mode for the case C = i and = 0.2. 



arises from the physical process 7 —>■ e^e . In the world-hne instanton interpretation 
this is described tunnehng process. 



Now that we have explained the origin of the singularities, we can continue with 
the analysis of ( [4.14| ). In order to produce a convergent integral, the Wick rotation 
T — > —iT can be performed as illustrated in Figure (H). Once Wick rotated, the 



X X X X — >► 



T 



Figure 8: Wick rotating the integration contour to the negative imaginary axis. The crosses 
represent branch point singularities. 



resulting integral has the form of an inverse Borel transform: 

,2 roo+ie jrp ^11 ^2 



ni{uj) = 1 - 



am 



dT 



n 



sinh Pi cosh A y sinh (3i cosh A 

(5.10) 

When CJ2 is imaginary, i.e. the Type II case, there are branch point singularities on 
the integration contour and one has to be careful to take the contour to lie above 
the real axis as indicated by the ie prescription. 



Since the T integral is of the form dT e ^f{T) where the function /(T) 

satisfies the reality condition /(T*)* = f{T), it follows that the imaginary part of 
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the integral is equal to | J^dT e"'^ f{T), where C is a contour that comes in from 
oo below the cut, goes round the first branch point singularity and then goes to oo 
above the cut, as illustrated in Figure (H). Hence, 



>: 

>; 



T 



Figure 9: After Wick rotation, the contour 6 that computes the imaginary part of ^^((jj). 



Im ni{u) 



am 



sinh Pi cosh Pi 



n 

1=1 



Pf 



sinh Pi cosh Pi 



(5.11) 



The imaginary part of the refractive index has an interesting interpretation be- 
cause it computes the probability for pair creation, 7 e~^e~. In fact the total 
cross-section per unit volume, or inverse mean free path, is 



^m.f.p. 



uj Imn(uj) 



(5.12) 



Notice that it is only non-vanishing in the Type II case when there are singularities 
on the real axis in the (Euclidean) T plane. Earlier we pointed out that these sin- 
gularities correspond to non-trivial classical loops with imaginary affine parameter. 
We can now identify these solutions as world-line instantons that describe the tun- 
nelling process 7 e~^e~. The fact that they occur only when (T2 is imaginary is 
natural. Remember, when (T2 is imaginary the null congruence is defocusing in the 
y"^ direction. This suggests the following intuitive picture: after a virtual e~^e~ pair 
is produced by tunnelling, the pair then follow diverging geodesies and become real 
particles. 



At low frequencies, the probability is dominated by the position of the first 
singularity at 



in Euclidean space, corresponding to the fundamental world-line instanton which 
looks exactly like the red loop in Figure (^ in Euclidean time. The Euclidean action 
of a zero mode is simply 5*^ = T, and for small u we can use the steepest decent 
method to approximate the ^ integral. The saddle-point is at ,^ = and hence 
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the leading order behaviour of ( p.ll ) for small uj will be of the form of an essential 
singularity: 

Lmuiiuj) ~ exp j — - . v^-^"^) 

u; 0-2 



5.2 The low frequency regime 



Low frequency means that uj'^R/rn^ ^1. As a consequence of this, the length of 
the loop ( p.l6 ), L ~ oj/m?^ is much smaller than the curvature scale: LR^/^ = 



LoR^^"^ /m? <^ 1. The leading order term in this limit will consequently be insensitive 
to the u dependence of hij{u) and so our result for this term is valid for all background 
metrics and not just ones which yield symmetric plane waves in the Penrose limit. 



To calculate the expansion in uj, we expand the Wick rotated integrand in ( p.lU| ) 
in powers of uo. The first term in the expansion is uo independent: 

and so there is no dispersion in this limit. Using (|3.3| ), this can be written in terms 
of the curvature, and the Newman- Penrose scalars, as 

/ N -, CI Ruu O! Ruiiu , /r^i 2\ 

1207r ra^ OOtt j^g-j 
JoOvr ^ 

for i = 1,2. In principle, this low frequency expression should follow from the terms 
in the low energy one-loop effective action of scalar QED that are quadratic in the 
field strength F and linear in the curvature. These terms have been calculated in 
spinor QED [1, 17, 18] but not, to our knowledge, in scalar QED.^° 

It is interesting to consider the higher terms in the frequency expansion in cer- 
tain particular examples. For example, for the case of a Type I conformally fiat 
background, ai = a2 = R^^"^, the velocity eigenstates for both polarizations have 



aK ^1 _ 71 u^R 428 fuj^Ry _ 15688 /u^Ry 

^'^^^ ~ ~ 27rm2 [l8 ~ 14175 ^ 189189 ~ 6891885 vT^J ^" 

(5.17) 

This series is divergent but alternating and this is correlated with the fact that it 
is Borel summable, with the sum being defined by the convergent integral in ( ^.10| ) 



^'^Even given the effective action, one must be very careful in simplifying with integration by 
parts because the on-shell photon wavefunction does not fall off at infinity. 
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which has no singularities on the real axis. Notice that ni{uj) is real to all orders 
in the expansion and since there are no cuts on the real axis the imaginary part 
vanishes, as is evident in ( |5.11| ). 



For the Type II Ricci flat background, cti 
superluminal at low frequencies with 



i(T2 = -R^/^, one polarization is 



ni{uj) = 1 — 



aR 



37 u'^R 34 /u;2^\2 



L45 28350 85995 



For the second, subluminal, polarization eigenstate. 



n2{oj) = 1 + 



aR 



1 



37 uj'^R 



34 



45 28350 m4 85995 V 



/u'^R\^ 



43 



135135 



/UJ^R\3 



43 



135135 V m 



(5.18) 



(5.19) 



The first series ( |5.18| ) is just the alternating version of ( |5.19| ). In both cases the Borel 
transforms have branch point singularities on the real axis and this is indicative of 



an imaginary part (5.14) which vanishes to all orders in the uj R/m expansion. 



5.3 The high frequency regime 

In the high-frequency limit uj'^R/mf" ^ 1, by re-scaling T m?T / {uj^{l — ^)) and 
expanding exp — m^T/(co'^(l — ^)) = l + -- -,we can show that the ni[ijj) approach 1 
like l/oj: 

n.M = l--^ + 0(!^) (5.20) 

where Cj is the integral 



oo+ie 



C=l 1 ("^^)-^ (521) 

* Jo sinh a^T cosh cTiT Aiy sinhVjT cosh a;T J ■ ^' ' 

Notice that the behaviour of the subleading term is softer than l/tu^. For the con- 
formally flat case with ai = a2 = R^^"^, the integral ( p.21| ) can be evaluated exactly 
by contour integration yielding 

for both 2 = 1,2. 

For the Ricci fiat Type II case (the vacuum gravitational wave), ai = ia2 = R^^^, 
although we cannot evaluate Cj analytically, there is an interesting relation 

C2 = -tCl , (5.23) 
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that follows from the definition of the integrals. A numerical evaluation in this case 
gives 

Ci = (0.22 - 0.014i)i?i/2 ^ ^ ^0.014 - 0.22i)R^/'^ , (5.24) 

which implies that both polarization states are superluminal at high frequencies. 
Hence, n2{uj) must change from being greater than 1 to less than 1 at some interme- 
diate frequency. 



5.4 Numerical analysis 



n-l 




Figure 10: The behaviour of ni(w) — 1 = n2{uj) — 1, in units of aR/{2TTm'^), as a function of 
^logiij^ R/m^ for the Type I conformaUy flat case ai — a2 = R^^^. The intercept ni{0) — 1 = 
~ -0.056. 



Type I: In this case, the integrand ( |5.10|) is regular on the real axis and so the 
resulting refractive indices are real and there is no pair creation. Figure (|T0|) shows a 
numerical evaluation n^uj) for the conformally flat background with ai = a2 = R^^^. 
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Figure 11: (a) The behaviour of Reni{uj) — 1 (i = 1 red, i = 2 green), in units of ai?/(27rm^), 
as a function of ^logcu^R/m'^ for the Type II Ricci flat case (vacuum gravitational wave) 
CTi = 1(72 = R^^^. Notice that the intercepts Reni(O) and Ren2(0) he the equal amount ^ ~ 0.023 
below and above 1, respectively, in accordance with the polarization sum rule [31]. The resolu- 
tion is not sufficient to show that the low-frequency subluminal photon becomes superluminal at 
high frequency, (b) A close-up of the region where Ren2(cj) — 1 changes sign signalling that the 
subluminal photon becomes superluminal at sufficiently high frequency. 
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Figure 12: The behaviour of (a) Inini(aj) and (b) Imn2(ti') for the superluminal polarization 
state, in units of aR/ (27rm^), as a function of i logijj^i?/m'* for the Type II Ricci fiat case (vacuum 
gravitational wave) cri — ia2 — R^^^. Notice that the subluminal polarization state, which is aligned 
with the defocusing direction in the null congruence, has a much larger value of Imn(cj) than the 
superluminal state. 



Type 11: In this case, the integrand ( |5.10| ) has branch point singularities on the 
real axis. From the point of view of a numerical evaluation, it is therefore not 
useful to perform the Wick rotation. A useful alternative is to perform a "half" 
Wick rotation by rotating the contour of ( [4.14| ) to lie along T ^ (1 — i)T/^/2. The 
resulting integral is convergent and can then be evaluated numerically. We find that 
the refractive indices have both a real and imaginary part, as we anticipated earlier. 
Figures (0) and (|l^) show the real and imaginary parts of ni{uj) for the example of a 
Ricci flat background with ai = ia2 = R^^^. Notice that the subluminal polarization 
state 71,2 (cu) behaves superficially like our simple model of a dissipative medium in 
Figure (B|). 



6. Micro- Causality and the Kramers-Kronig Relation 



Before we analyse our curved spacetime result, let us first consider the simple model 
of a dissipative medium. In that case, from ( |5.2| ) we see that n(to') has simple poles 
in the lower-half plane at uj = ituo — 1^/2 (for 7 <C cuq). Hence, n{uj) is analytic in 
the upper-half plane and the Kramers-Kronig relation is trivially satisfied. To see 
this, consider duj/uon{uj) for a contour along the real axis, jumping over the simple 
pole at u; = 0, completed by the large semi-circle in the upper-half plane, illustrated 
in Figure (|T3p. If n{ijj) is analytic in the upper-half plane, the total integral is and 
so: 



= / — n{uj) - mn{0) + V — n{uj) 

^semi-circle ^ J —oo ^ 

= 7ri(n(oo) — n(0)) + P f — n{uj) . 

J-00 ^ 



(6.1) 
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Taking the imaginary part, and assuming that n{oo) and n{0) are real and that 
Imn(a;) is an odd function, immediately yields (|1 . 1|) . It is a simple matter to check 
the relation explicitly for the dissipative model (|5.2|). 




Figure 13: The integration contour for § du) n{uj) / u) used in the proof of the KK relation for (a) 
the simple dissipative model with poles lying under the real axis (b) conformally flat case with poles 
on the imaginary axis. 



For curved spacetime this argument fails because there are singularities on the 
imaginary axis which have to be included in ( |6.1|) , as illustrated in Figure. ([T3|). For 
example, for the conformally flat Type I case, cxi = (72 = i?^''^, the singularities are 
poles whose residues must be included: 

m{n{oo) — n(0)) +vf —n{uj) = pole contribution . (6.2) 

Since in this case lm.n{uj) = 0, and including the contribution from the poles on the 
imaginary axis, (|6.2|) becomes 

Ren(O) -Ren(oo) = dT e-^ {^{l - S)f S^Res f{iTin/2) , (6.3) 

where we have defined the function 

f{x) = (1 — x^/(sinh^ x cosh^ a;))/x^ . (6-4) 

The residue sum can be regularized by considering f{x)e^"'^ and taking a — > at the 
end. The result is 

aR 

Ren(O) -Ren(oo) = -— (6.5) 

367rm^ 

which is in perfect agreement with ( |5.17| ) and ( |5.20| ). Notice that we have established 
this result by interchanging the order of the uj and T integrals in which case the 



^^For more general examples, the singularities are branch points and the integration contour has 
to come around them from e + too down the imaginary axis before going back to — e + too and 
completing the semi-circle. 
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singularities appear as poles on the imaginary axis. However, if we perform the T 
integral first, then the singularities become a branch cut in uj from to oo in the 
upper half plane. 

The fact that n{uj) is not analytic in the upper-half plane is intimately connected 
with the issue of micro-causality, as we now explain. In our simple model of a 
dissipative medium, the Fourier transform of the susceptibility, x{'^) — ~ 
l)/(47r), 

/oo 
du;e--*x(^) , (6.6) 
-oo 

plays the role of a response function: D{t) = E(t) + J dt' G{t — t')E{t'). In the 
simple model, n{uj) and hence x(c^) is analytic in the upper-half plane and so, when 
t < 0, we can compute the u integral by completing the contour with a semi-circle at 
infinity in the upper-half plane. Since there are no singularities, the integral vanishes 
implying G{t) — 0: cause precedes effect. Taking the explicit Fourier transform, we 
have 

G{t) = ^e-^'/^ sm{ujot)e{t) . (6.7) 

In the curved spacetime case, n{uj) is not analytic in the upper-half plane and so it 
implies that the analogue of G{t) will be non-vanishing for t < 0. 

We now place this simple analysis in the context of relativistic QFT, where 
response functions are more properly understand in terms of (retarded) propagators. 
The one-loop vacuum polarization n^"'°°P contributes to the propagator via A = 
^tree _ ^tree-Qi-ioop^tree _|_ . . . ^ real spacc picturc, the issue of micro-causality 
rests on the fact that the retarded propagator Ai.ct(ic) is only non- vanishing in, or on, 
the forward light cone.^^ In prosaic language, an external source can only infiuence 
the fields in the future. For instance, in the present context the real space tree-level 
retarded propagator A*g^^ is only non- vanishing on the forward light cone. However, 
what about the one- loop correction? Notice that we have only calculated H^"'°°f'(cL)) 
on-shell in momentum space and this means that we do not have access to the 
complete one-loop real space propagator. However, we can perform the Fourier 
transform with respect to cu, which determines the propagator as a function of the 
null coordinate v: 

/oo 
(ia;e-''"^H^-^°°P(u;) . (6.8) 
-oo 

This is a retarded quantity if the integration contour is taken to avoid singularities 
by veering into the upper-half plane, when v < 0, and the lower half plane, when 
V > 0. For QFT in flat spacetime, H^"'°°P(a;) is analytic in the upper-half plane 

^^When we talk in the following about the "light cone" we mean the geometrical null surface 
defined by the metric gr^,^. 
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and so when v < one computes the u integral by completing the contour with a 
semi-circle at infinity in the upper-half plane. Since there are no singularities in the 
upper-half plane, the integral vanishes and consequently njg^°°'^(f) = for w < 0. 
This is consistent with the fact that the region f < lies outside the forward light 
cone. Hence, in this case micro-causality is preserved as a consequence of analyticity 
in the upper-half plane in frequency space. In curved spacetime, on the contrary, 
n^"^°°P(ci;) is not analytic in the upper-half plane and consequently it would seem 
that the one-loop retarded propagator Hl'J^'^'^^v) must receive contributions from the 



region v < which lies outside the forward light cone. See Figure (|Tj) 




Figure 14: Including vacuum polarization effects, the photon momentum k may lie outside the 
forward light cone {u > 0, v > 0) of its original null geodesic v — 0. The potential violation of 
micro-causality implies that the retarded propagator is non-vanishing even for v < (the shaded 
area), which lies outside the forward light cone. 



The idea here, is that uo is identified with one of the lightcone momenta When 
the photon is on-shell at tree level, the other component vanishes, p_ = 0. Now by 
giving p_ a small positive imaginary part a non-analyticity in the upper-half p+ plane 
would be in the region Imp+/Imp_ > which means that the retarded propagator 
is non-vanishing outside the lightcone. However the loophole in this reasoning is 
that once the T and ^ integrals have been performed, the non-analyticities of n{uj) 
as a function oi u = arises as branch cuts joining p+ = to j9+ = ioc. When 
taken off-shell, it may be that the branch cut slips into the causally safe region 
Imp+/Imp_ < 0. The only way to really settle this issue is to perform a calculation 
of the vacuum polarization with the momentum off-shell [39].^^ 

One final point to emphasize is that in the Type II examples where the imaginary 
part of the refractive index is non- vanishing, it is positive as expected from the optical 

"^^Note added: We have now completed a full off-shell calculation of the vacuum polarization 
tensor and find that precisely this behaviour occurs. The branch point at = ioo is shifted to 
— rr? Ip_ with Imp_|_/Imp_ < 0. Full details will be presented elsewhere [40]. 
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theorem^^ which relates it to the total cross-section per unit volume, or inverse 
mean free path as in (|5.12| ). In this sense, spacetime acts as an ordinary dissipative 
medium. However, unlike the simple dissipative model where the imaginary part 
of the refractive index falls off as l/cu^, for curved spacetime the imaginary part 
falls off as OiB}!"^ jijj. This implies that the mean free path saturates to a constant 



7. Remarks on General Backgrounds 



In this section, we consider some of the features of our analysis which apply in more 
general non-symmetric plane wave backgrounds when hij{u) = Ruiju{u) does not 
take the special form ( |3.1D . For the moment, we shall assume that the Penrose 
limit is non-trivial, i.e. not flat space. Our purpose is to highlight how the resulting 
behaviour of the refractive index depends to a large extent on the properties of the 
null congruence. 

At the beginning of Section 5, we found that the analytic structure of the inte- 
grand ( |4.14|) - more precisely the positions of the singularities - could be traced to 



the existence of zero modes of the fluctuation equations. In the general case, the 
equation for these zero modes is 

f + iL%,y^ + '^n.^y^ir - - '^^i.y^Sir) = -C . (7.1) 

Considering the solutions of these 2nd order equations in the two regions < r < ^ 
and ^ < T < 1, and including C and T, there are generically nine unknowns to be 
fixed, up to overall scaling of the solution. At r = and r = ^ there are a total of 
eight boundary conditions on ?/* and y*, with the constraint rfry* = providing a 
ninth condition. In general, we therefore expect solutions to exist only for particular 
values of T. This is exactly what we found for the symmetric case where the special 
values of T are given in (|5.4|) . Although, for n even we found two zero modes rather 
than the one expected. 

As we have seen, the existence of the zero modes is related to the behaviour of 
the null congruence. This can be made more concrete by looking at the special case 
when ^ = I (which is, in any case, picked out by the saddle-point method described 
above). In this case, by symmetry we expect a solution of the form that we found 
in the symmetric plane wave case, ( |5.5| ), but where the transverse geodesic deviation 



^^It is amusing to note that here we are applying the optical theorem in its original context of 
the refractive index from which the name of the theorem is derived. 
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vector y''{u) satisfies the equation for a Jacobi field along the geodesic 7: 



subject to the boundary conditions ?/*(±Mo) = 0. Because of the latter, the source 
terms vanish and so the derivatives dy^/dr must be continuous at r = 0(= 1) and 
r = ^, a fact that follows directly from the ansatz ( ^.5|) . In addition, the constraint 
dry^ = is automatically satisfied and the Lagrange multiplier C vanishes. The 
classical solution is then u = u{t), (with a slight abuse of notation) ?/*(r) = y^uij)) 
and v{t) satisfies its own geodesic equation. Since ( |7.2|) is a second order linear equa- 
tion there will in general be solutions only for particular values of Uq = ujT/{Am^). 
Since y'^{u) vanishes at m = i-uo, at least when the special values of T are real, these 
points are precisely conjugate points along 7. Hence, the existence of zero modes 
(for real values of T) is tied directly to the existence of conjugate points. Notice, 
however, that the special values of T for which zero modes exist are not necessarily 
real. This is exactly what happens in the Type II examples, where the singularities 
corresponding to world-line instantons have imaginary T. 

The zero modes dictate the analytic structure of the T integral which in turn 
determines the nature of the physics. In particular, the singularities along the real T 
axis play a prominent role because, as we have seen, they are responsible for the non- 
trivial analytic structure of n{uj). Moreover, as we have argued above, zero modes 
for real T correspond directly to the existence of conjugate points along 7. But the 
existence of these conjugate points, as explained in Section 3, is generic. Therefore 
we are led to the following conclusion: 

Conclusion: violations of analyticity and the Kramers- Kronig relation 
are generic and can be traced to the focusing nature of null geodesies and 
the existence of conjugate points implied by the null energy condition. 

As already mentioned, the Penrose limit is ideally suited to the analysis of photon 
propagation in arbitrary background spacetimes. Many of the characteristic features 
of superluminal low-frequency propagation previously found in specific examples, 
including Schwarzschild, Reissner- Nordstrom and Kerr black holes [1,15,31] as well 
as gravitational waves [1, 16], can be seen directly in the Penrose limit. For example, 
a maximally symmetric spacetime such as de Sitter has vanishing $00 aiid ^'o and the 
low- frequency phase velocity Wph(O) receives no correction from vacuum polarization. 
Using our formalism, we see immediately that at leading order in Rjvr? this result 
holds for all frequencies since the Penrose limit of a maximally symmetric spacetime 
is fiat [23]. 



- 35 - 



In Schwarzschild spacetime, we have previously found that while a photon follow- 
ing a general null geodesic may experience a superluminal shift in fph(O), the effect 
vanishes for purely radial geodesies. (In fact, this remains true for photons following 
principal null geodesies [30] for any Petrov type D spacetime such as Schwarzschild 
or Kerr, again due to the vanishing of the corresponding $00 and \E'o.) This is clear 
in our formalism. The Penrose plane wave limit for the Schwarzschild metric is, in 
Brinkmann coordinates, 

ds' = 2dudv+^-^{{yr - {y'f)du' - {dy^f - {dy')\ (7.3) 

where L specifies the angular momentum and r{u) is given by the solution of the 
geodesic equation. We see immediately that for radial trajectories the Penrose limit 
is fiat and so, at least at 0{R/rn?), the phase velocity fph(i^) remains equal to c 
for all frequencies, not just in the low-frequency limit. Clearly, in such cases where 
the Penrose limit is fiat, the expansion (|2.30| ) gives a systematic way to go beyond 
leading order in R/m?. An interesting feature is the existence of a "peeling theo- 
rem" [23] , whereby successive orders in the Penrose expansion involve the curvatures 

This gives a first glance at the power of the Penrose plane wave geometry com- 
bined with the world-line sigma model approach. Moreover, other general features 
of null congruences will play an important role. For example, we have been implic- 
itly assuming that the geodesies are complete so that the affine parameter varies 
from —00 to +00. However, there are spacetimes where certain null geodesies are 
incomplete and the affine parameter has a finite limiting value. This usually signals 
the existence of a spacetime singularity, as for example in the case of Schwarzschild 
orbits for L less than a critical value, where the Penrose limit becomes singular [23]. 
Clearly, this can affect the zero modes in the sigma model and therefore the singu- 
larities and asymptotic behaviour of the refractive index. The role of horizons in 
relation to the Penrose limit also deserves investigation. All of these issues will be 
considered in detail elsewhere. 



8. Conclusions 

In this paper, we have for the first time evaluated the non-perturbative frequency 
dependence of the vacuum polarization for QED in curved spacetime and determined 
the corresponding refractive index for photon propagation. In so doing, we have 
resolved the outstanding problem in "quantum gravitational optics" [3,4], viz. how 
to reconcile the prediction of a superluminal phase velocity at low frequency with 
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causality. Remarkably, the resolution involves the violation of analyticity calling into 
question micro- causality in curved spacetime. 

These results have been achieved by combining two powerful techniques: (i) the 
world-line sigma model, which enables the non-perturbative frequency dependence 
of the vacuum polarization to be evaluated by a saddle-point expansion around a 
geometrically motivated classical solution, and (ii) the Penrose plane wave limit, 
which encodes the relevant tidal effects of spacetime in the neighbourhood of the 
original null geodesic traced by the photon. 

The form of the refractive index reflects the nature of the background spacetime. 
We identify two classes. In Type I backgrounds, which include conformally flat 
spacetimes,^^ both photon polarizations are superluminal at low frequencies, but the 
phase velocity approaches c at high frequency. The imaginary part of the refractive 
index vanishes. In Type II backgrounds, which include Ricci flat spacetimes, photon 
propagation may display birefringence with one superluminal and one subluminal 
polarization at low frequency. In both cases, however, the high frequency phase 
velocity is c. The refractive index develops an imaginary part, indicating a non- 
zero probability for pair creation, 7 e+e~. Since the high-frequency limit of the 
phase velocity is identified with the wavefront velocity f^f, which is the "speed of 
light" relevant for causality, we see explicitly how superluminal propagation in the 
low-frequency theory is compatible with causality. 

Although these results were obtained using the Penrose limit in locally sym- 
metric spacetimes, they are expected to be generally true. The reason is that the 
analytic properties of the refractive index can be related in the world-line sigma 
model formalism to general results in the theory of null congruences. In particular, 
the distinction between Type I and Type II spacetimes is whether the null geodesies 
in the congruence focus in both transverse directions (Type I), or focus in one and 
defocus in the other (Type II). The result that at least one direction is focusing is 
a consequence of the null energy condition. The presence of a focusing direction in 
the congruence then implies the existence of conjugate points, which leads to the 
existence of zero modes and ultimately yields poles in the refractive index in the 
upper-half complex plane, violating the analyticity assumptions used to derive the 
Kramers-Kronig dispersion relation. The violation of this dispersion relation in turn 
allows n(oo) > n(0) and removes the apparent paradox of having a superluminal 
phase velocity fph(O) > c while the wavefront velocity Wwf = t'ph(oo) = c. 

This is potentially the most far-reaching conclusion of this paper. The null 

^^Note that the Penrose Umit of a conformally flat spacetime is also conformally flat. Similarly 
for Ricci flat, and also locally symmetric, spacetimes [23]. 
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energy condition and the general relativistic theory of null congruences necessarily 
imply a non-analyticity of the refractive index, although the full implications of this 
for micro- causality and the other axioms of S-matrix theory will only follow from an 
off-shell extension of the calculation. 

The loss of analyticity in or more generally in forward scattering ampli- 

tudes, also has important implications for the idea that constraints may be placed 

on the parameters of a low-energy effective field theory by the requirement that it 
admits a consistent UV completion [4, 37, 38] . These constraints are typically de- 
rived either by requiring the absence of superluminal effects in the low-energy theory 
or assuming analyticity in dispersion relations involving forward scattering ampli- 
tudes. While these remain valid in flat spacetime, we have shown that they are not 
applicable to fundamental UV theories involving gravity, including string theory. 

The full imphcations of the calculation of the refractive index and the issues of 
causality and micro-causality remain to be explored, especially in relation to horizons 
and singularities. The significance of the UV-IR mixing whereby the high-frequency 
limit probes the global properties of the null geodesic congruence also deserves to 
be better understood. What is clear, however, is that the results described here will 
have a significant impact on our understanding of quantum field theories involving 
gravity. 

We would like to thank Asad Naqvi for many useful conversations and Sergei 
Dubovsky, Alberto Nicolis, Enrico Trincherini and Giovanni Villadoro for pointing 
out the necessity of working off-shell in order to completely settle the question of 
micro-causality. TJH would also like to thank Massimo Porrati for a helpful discus- 
sions and Fiorenzo BastianeUi for explaining some details of his work on the world-hne 
formalism. This work was supported in part by PPARC grant PP/D507407/1. 

Appendix A: Power Counting 

In this Appendix, we prove in an alternative way one of the key results of this 
paper: that each loop in the world-hne QFT comes with a power of Rju? and so 
loops are suppressed in the limit of weak curvature R <^ m^. In order to assess the 
behaviour of a given graph in perturbation theory, it is useful to re-scale T — > T jrv? 
and then r — > tT and x — > \/Tx so that the world-line action can be split as 




(A.l) 
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where a typical term in 5'pert arises from expanding the metric around flat space at 
the point xq = 0; schematically, 

1 / P \ "/^ 
dr[-,) (A.2) 

where i?"/^ denotes powers of the Riemann tensor and its derivatives of mass dimen- 
sions n. The vertex behaves as (i?/m^)"/^ and has n + 2 legs. In addition, we have 
the exponential factors uQ which we can view as additional vertices of the form 



Consider such a graph with E external legs, I internal legs and V vertices. If the 



graphs consists of A^^ vertices of the form ( |A.2|) and Sn vertices of the form (|A.3|) , 
then 

Y,{{n + 2)N^ + {n + l)Sr)=2I + E , ^ = ^ (iV„ + 5„) . (A.4) 

n n 

The graph behaves as 

^E„^?n^-E„(«JVn + (n+l)5„)^E„ri(JVn+5„)/2 ^ / j (~^) • (A-S) 

Now we use the topological identity, L = I — V + 1, where L is the number of loops, 
to equate this to 

So each loop brings a factor of R/m^. For example, the partition function Z has 
E = Q and since the tree-level contribution is the classical action for the saddle point 
which vanishes, the leading order term comes from one loop and is an arbitrary 
function of u'^R/mf'. The expansion around the classical saddle-point solution sums 
up all the one-loop graphs with arbitrary uj insertions. The leading order contribution 
to the Green's function piece, which has E = 2, comes from tree level. Once again, 
the expansion around the classical saddle-point solution sums up all these tree graphs 
with arbitrary u insertions. 
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